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We propose a new approach to extract the wave functions of resonances by the bound 
state approximation which gives the mixed states of the resonance components and the 
continuum ones. In our approach, on the basis of the method of analytic continuation in 
the coupling constant (ACCC), we construct Pade rational function by adopting the positive 
energies as well as the negative ones. We report the result of the application of this new 
method to the second 2 + state of 12 C which was studied with the ACCC method in our 
previous work. It is found that the resonance parameters obtained by the ACCC method 
are well reproduced by the new method. Some advantages over the ACCC method are also 
shown. 

§1. Introduction 

In order to study broad resonance states, it is well known that the bound state 
approximation does not work well. Among the methods which can treat resonance 
states, we recently see many studies which are based on the ACCC method^ 
and complex scaling methodP^^ In the case of the ACCC method, however, one can 
only calculate the energy and width of the resonance but to get the resonance wave 
function is almost impossible. On the other hand, in the case of the complex scaling 
method, one can in principle get the wave function in addition to the resonance 
energy and width, which is however not easy in practice. Moreover, the wave function 
which is to be obtained by the so-called back rotation applied to the complex-rotated 
wave function is the so-called Gamow state which is exponentially diverging in the 
outside spatial region and is not easy to handle. 

The purpose of this paper is to propose an extended version of bound state 
approximation to the description of resonances which provides us with a new method 
of the calculation of the energy, the width and the wave function of the resonance. 
In our extended version of the bound state approximation, we first perform the 
calculation of usual bound state approximation and obtain the energy eigen-states. 
Since the energy eigen-states consist of the mixture of the resonance wave function 
and the continuum state components, we then extract the resonance component by 
superposing energy eigen-states. For this procedure we need to know approximately 
how is the mixing between resonance states and continuum states, and this knowledge 
is obtained by using the pseudo potential just like in the ACCC approach. The 
difference of our method from the original ACCC method is that we can also keep 
the coupling constant of the pseudo potential so weak that the resonance state does 
not become the bound state, while in the original ACCC method it is inherently 
necessary to use the strong coupling constant which converts the resonance state 
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to the bound state. The first part of our method, namely the extraction of the 
resonance state component from energy eigen-states, was already utilized recently 
by the present authors in RefP Here in this paper, we give a full explanation of our 
method and the calculational method of the resonance width. 

The content of this paper is as follows. In £J21 we explain our new treatment 
of the resonance state in the bound state approximation with the aid of the pseudo 
potential. It consists of two parts. The first part is the extraction of the resonance 
wave function from the energy eigen-state wave functions obtained by the bound 
state approximation. The second part is the calculation of the resonance width. 
In |21 we show the numerical example of the application of this new method. The 
application is done to the second 2 + state of 12 C which was recently observed at 
2.6 ± 0.3 MeV above the 3a breakup threshold with r = 1.0 ± 0.3 MeVP Since 
this state was investigated by using the ACCC method in the previous paper by the 
present authors^ we here compare the results of the new ACCC and the previous 
one. We show that the complex energy function obtained by our method is in good 
agreement with the one of the original ACCC method in a wide coupling constant 
region, and therefore our method gives the energy and width which well agree with 
those obtained by the ACCC method. £0]is for conclusion. 

§2. New Treatment of Resonances with Bound State Approximation by 

the Aid of Pseudo Potential 

2.1. Extraction of resonance wave function 

The calculation of the resonance state in bound state approximation is usu- 
ally made by diagonalizing the Hamiltonian by using a finite number of basis wave 
functions which are square integrable. The calculated positive energy eigen-states 
which are linear combinations of basis wave functions are divided into resonance 
states and continuum states. There are many practical methods to make this divi- 
sion, namely the selection of the resonance states among positive energy eigen-states. 
One method is to introduce an attractive pseudo potential V which is added to the 
original Hamiltonian H, 

H'{5) = H + 5 xV, (2-1) 

where 5 is a coupling constant to vary the strength of the pseudo potential. We 
diagonalise this new Hamiltonian H'(5) by using the same set of basis wave functions 
as used for H. As we strengthen the coupling constant 5 from the physical value, 
5 = 0, the eigen-energy of any resonance state gets lower and finally the resonance 
state changes into a bound state. On the other hand continuum states show almost 
no change in their eigen-energies for the change of 5. 

These behaviors of the eigen-energies as functions of 5 are shown schematically 
in Fig. ^ In this figure, the curve with label (b) around 5 = shows the eigen-energy 
of a resonance state denoted as ^b(^) and goes down as 5 increases. This curve (b) 
then crosses with the curve (a) which is flat around 5 = and is the eigen-energy 
of a continuum state denoted as tpa.(5). After the closest approach of the two curves 
(b) and (a), the curve with label (b') which is the smooth continuation of the curve 
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Fig. 1. Eigen-energies of Hamiltonian (12-111 are schematically shown as functions of 8. There are 
three crossing points near which the mixing of continuum and resonance components occurs. 
See the text for the detailed explanation. 



(a) and hence is disconnected with the curve (b) represents the eigen-energy of the 
resonance state tpb(8). Thus, as usual, we say that the curves (b) and (a) have a 
crossing. The eigen-energy of the continuum state i/j a (5) which is first represented by 
the curve (a) continues to be almost flat except near crossing points with resonance 
curves and converts into the curve (a'). In the case of the resonance state V'b(^) 
at 5 = 0, we can regard that this state contains a negligible amount of mixture of 
continuum state component because this state is located far from the crossing points 
with continuum states. 

However, in the case of the resonance state ipd(5) on the resonance state curve 
(d) in Fig. ^ its crossing point with the continuum state curve (c) is very close 
to ipd($) at 5 = 0. Therefore ipd($) at (5 = should contain sizable mixture of the 
continuum state. This kind of situation is usual for a resonance state with broad 
width, and the resonance wave function obtained by the bound state approximation 
is not so much trustworthy because of sizable mixture of the continuum state. 

In order to resolve this difficulty, we note the utility of the continuum state wave 
functions obtained by the bound state approximation in the neighborhood of the 
resonance state, because they should contain sizable mixture of the resonance state 
component. In our present simple example of Fig. ^ we notice the utility of the 
continuum state wave function tp c (5). We can express both states, V'd(^) and ip c (5), 
as linear combinations of pure resonance state and pure continuum state &c($)', 
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MS) = a& d (6) + pP c (8), MS) = -PMS) + aV c (8). (2-2) 

The extraction method of the pure resonance component !?d(<5) which we propose is 
obtained by noticing that the density extension of a resonance state is far smaller than 
those of continuum states. Namely we can determine a and [3 from the requirement, 

A 

<!Pd(*)| 2>?|!Pd(<5)) = smallest, (2-3) 

MS) = »MS) ~ PMS), (24) 

where A is the mass number. This requirement is satisfied by diagonalizing the 
operator J2j=i r j by the basis states, ipd(5) and tp c (S), and by choosing (a,— (3) to 
be the eigen-vector belonging to the smallest eigen-value. 

When the resonance component is distributed in several energy eigen-states, we 
need to diagonalize the operator Ylf=i r j by such several energy eigen-states. When 
the energy eigen-states contain two resonance components, the diagonalization of the 
operator Ylf=i r ] & ves us two small eigen-values. Two resonance states are obtained 
by diagonalizing the Hamiltonian H'{8) with two wave functions which give the two 
smallest eigen-values of the operator Ylf=i r ]- 

We are now able to calculate approximate values of resonance energies along the 
resonance energy curve. If a point on the resonance state curve is out of the influence 
of any crossing points with continuum state curves, we can regard that the energy of 
the point is just the resonance energy. On the other hand, if a point on the resonance 
state curve is located within the influence of some near-by crossing point, we extract 
resonance wave function by the procedure explained in the previous subsection and 
then calculate the expectation value of the Hamiltonian H'{8) with the extracted 
wave function which is just the resonance energy at the point. 

2.2. Beyond the ACCC method 

In the ACCC method, the complex eigen-energy E{5) of the Hamiltonian H'{§) 
is approximated by Pade rational function, 

9 N 

_ . C + C\X + c 2 x z -I 1- cnx" 

[N ' M] ~ \ + d 1 x + d 2X i + --- + d M x M 

E(5) - E th (5) = kf NM] , 

where E t h(S) is the threshold energy, x = y/5 — 5q, and 5q is determined by the 
reqirement E(5q) = E t h(So). cq is reserved to take care of the possible errors in 
determining 5q in practical calculations, though it has to be zero in principle. The 
(N + M) coefficients, c\ ~ cat, d\ ~ du, are determined by calculating E{8) — 
-E t h(<5) for (N + M) different values of 5 in the bound state region, 5 > Sq, namely 
the region where E(S) — E t h(5) < 0. 

What we want to know is the complex eigen-energy E{5) at 5 = 0. However, in 
the ACCC method we have to use the information in the bound state region where 
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5 is larger than 5$ and is far from 5 = 0. It is the reason why we need high accuracy 
in the calculation with ACCC. 

As we explained in the previous subsection we are now able to calculate ap- 
proximate values of resonance energies along the resonance state curve. The knowl- 
edge of the resonance energies enables us to determine the (N + M) coefficients, 
ci ~ cat, d\ rsj oIm, as we explain below. Once the Pade rational function is obtained, 
we can calculate the resonance widths by using it. The merit of our method is that 
we can construct the Pade rational function not only by using the bound state in- 
formation with 5 > 5q but by using the information in the resonance state region 
where 5 is smaller than 5o and is near 5 = 0. 

The explicit procedure to determine the (N + M) coefficients, c\ ~ cjy,di ~ oIm, 
is as follows. We calculate the resonance energies 1ZeE(5) for L different values of 5, 
5 m ,m = 1 ~ L. The constant L is an arbitrary integer. Then we have the following 
relation, 

2(neE(5 m ) - E th (5 m )) = (k [N>M] ) 2 + (k* [NM] ) 2 

X m = \/5m - 5q (5 m > 5 ), 
X m = ~W 5 - 5 m (5m < 5 )- 



(2-7) 
(2-8) 



The determination of the (N + M) coefficients, c\ ~ CN,d\ ~ dm, by solving 
these simultaneous equations can be made by using the following frictional cool- 
ing methocP^ used widely in the works with AMD (antisymmetrized molecular 
dynamics) theory JSBE) 



dt° m 



!/'■ I 



dF 

dc m . 
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F=J2 [ 2 (KeE(5 m ) - E th (5 m )) + (- 
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Sj=0 c j( x m,y 

1 + T,jLidj(x m y 



(2-9) 



(2-10) 



Here the (N + M) coefficients, c\ ~ cjy,di ~ dM, are regarded as being functions 
of time t and /x is an arbitrary non-zero real number. It is easy to show that F 
decreases as time develops, (d/dt)F < 0. Hence after sufficient time steps, we obtain 
the minimum value of F, which is zero. Needless to say, F = is equivalent to that 
the simultaneous equations of Eq. (|2-7[) are satisfied. It is to be noted that the initial 
values of ci ~ cat, di ~ dM for solving the frictional cooling equation of Ea. (|2-9(l are 
arbitrary. 
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§3. Application to the state of 12 C 

In the previous paper by present authors,® the 2^ state of 12 C was studied with 
the following type of wave function, 



v 3a 



03) = Pj= 2 ^[nexp { - ( ^{X ix - X Gx f + -^(X iy - X Gy ) 2 
i=i x v 

+^(X iz -X Gz f)}<j ) {a i )\. (3-1) 

Here, B 2 = b 2 + 2(3 2 , (k = x,y,z). 0(oi) oc exp[-(l/26 2 ) £j(r iW - X;) 2 ] with 
Xi = (1/4) £ • r j(i) an d = 4(i — 1) + j is the internal wave function of the i-th 
alpha cluster, Xq is the total center-of-mass coordinate and Pj=2 is the projection 
operator of the angular momentum onto J = 2. The wave function of Eq. (|3-l|) 

is just the 3a case of our general na Bose-condensed wave function proposed re- 
centlyPJCm The 

eigen-energies in bound state approximation are obtained by 
solving the following Hill- Wheeler equation, 

£ ($ J 3 : 2 ((3)\(H'(5) - E x (5))\$i: 2 ((3'))f J x = 2 ((3') = 0. (3-2) 
0' 

The solution of the above equation is given as, 

^/ =2 W = ^//= 2 (/3,^= 2 (/3). (3-3) 
P 

The microscopic Hamiltonian is expressed as 

H = T - T G + V N + V C . (3-4) 

Here T stands for the kinetic energy, Tq the total center-of-mass kinetic energy, Vq 
the Coulomb energy between protons, and Vn the effective two-nucleon force. In 
this paper, the Volkov No.2 forced with Majorana exchange mixture M = 0.59 is 
adopted as Vn and the following form as the pseudo potential V in H'(5) = H+8xV, 

ij 



r—<ww\^-^£-.). (35) 

We should note that the same Hamiltonian H'(5) as in RefP is adopted in this 
paper, so that we can check whether our method works well or not by comparing 
the results of the new ACCC with those of the original ACCC used in Ref)^ 

In Fig. I2 the complex energy function k 2 7 « (5) for N = M = 7 calculated 
by using the extended version of ACCC method is drawn, together with that of 
the original ACCC method which is denoted as k\ c {5). The k\n{5) was carefully 
determined by using the Pade approximant of N = M = 8 in RefP* Dotted curves 
denote the eigen-energy levels obtained by solving the Hill- Wheeler equation Eq. 
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Fig. 2. The comparison between the complex energy functions which are obtained by the new ACCC 
method and the original one. These functions are denoted by solid curves. The difference 
between the real parts of the complex energy functions obtained by both methods is almost 
invisible. Dotted curves show the energy-eigen values which are the solutions of Hill- Wheeler 
equation Eq. (13-21 . Cross points denote the energy levels obtained by the extraction method 
explained in fl2.ll where the energy-eigen values denoted by open circles in 5 < 0.004 are used. 
The energy levels denoted by open circles in 5 > 0.004 and cross points are adopted to construct 
the complex energy function expressed by the Pade approximant. 



(|H-2|) . The eigen-energy levels which are used in order to determine the coefficients, 
ci ~ C7, d\ ~ d-j of the Pade rational function Eq. l|2-5j) are denoted by open 
circles. We superpose the denoted levels by open circles in order to derive the correct 
resonance states by the extended version of bound state approximation explained in 
^2.11 which we simply call the extraction method in the following. This is done at 
each position of 5 which is less than 0.004 in the resonance region. The obtained 
energy levels by the extraction method are denoted by crosses. 

We see that the energy levels which are given by the extraction method and 
denoted by crosses are almost exactly on the trajectory 7£e(£r^ c (<5)). This means 
that the energy levels given by the extraction method well reproduce the resonance 
energies predicted by the original ACCC method in the wide region of 5. It is to 
be noted that since in this 6 region less than 0.004 the influence of the near-by 
crossing point is not negligible, without the extraction method the eigen-energy of 
usual bound state approximation cannot be on the trajectory TZe(kj iC (5)) . Thus we 
can say that our extraction method works effectively and therefore, the obtained 
wave function is more or less reasonable. We have made a detailed analysis of thus 
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obtained wave function of the state in RefP 

Reflecting the fact that all of the energy levels in the new ACCC method, i.e. 
the levels denoted by crosses for 5 < 0.004 and by open circles for 5 > 0.004, are on 
the trajectory lZe{k\ G {5)) the real part of the complex energy function TZe(k^ 7 « (5)) 

coincides with that of the original ACCC method TZe(k\ c (5)). The coincidence of 
the imaginary part Im(k^ 7 7 j(5)) with Tm(k\ c (5)) is not so good as the real part. 
Nevertheless, we can say the agreement of both trajectories is nice by taking into 
account the fact that the new and original ACCC method are usually followed by 
numerical errors to some extent. The calculated resonance energy TZe(k^ 7 ^(5 = 0)) = 

1.59 MeV, which is mesured from the 3a threshold, and the width —2Im(k? 77 JS = 
0)) = 0.15 MeV, while lZe{k\ c {5 = 0)) = 1.55 MeV and -21m{k\ c {8 = 0)) = 0.10 
MeV in the original ACCC method. We have checked that the other choices of 
[N, M] of the Pade rational function, namely N = M = 6, 8 resulted in almost 
the same trajectories for both real and imaginary parts. Incidentally, the observed 
binding energy and width for the 2^ state in 12 C are 2.6 db 0.3 MeV measured from 
the 3a threshold and 1.0 ± 0.3 MeV, respectively. It may be said that the theoretical 
values of the width is much smaller than the experimental value. We discussed 
this point in the previous work,® where it was shown that this occurs because the 
theoretical binding energy is lower than the experimental one and this makes the 
decay penetrability smaller. See RefP for further detailed discussion. 

We examined other input eigen-energy values or other number of the adopted 
eigen-energy levels in determining the coefficients of the Pade approximant. Except 
for the case that the eigen-energy values in bound state region are not used, the 
resultant complex energy function is not ill-behaved but qualitatively unique. The 
merit of the new ACCC method is that arbitrarily large number of eigen-energies can 
be used to determine the (N + M) coefficients of Pade rational function so far as the 
valid numerical accuracy is assured in the frictional cooling method. This should be 
compared with the original ACCC method which allows just (JV + JW) eigen-energies 
to be used. Owing to this advantage, we can obtain the correct trajectory of the real 
part of complex energy function with almost no ambiguity in the wide range of 6 up 
to 5 = 0, and therefore the imaginary part can necessarily result in almost unique 
trajectory at least in the adopted 5 region. It is well known that in the original 
ACCC method the divergence or ill-behavior easily occur in the obtained energy 
function unless sufficiently high numerical accuracy is assured. 

As mentioned above, in the new ACCC method we can construct the Pade 
rational function not only by using the bound state information with 5 > 5o but also 
by using the information in the resonance state region where 5 is smaller than 6q and 
is near 5 = 0. This allows us to determine the coefficients of Pade rational function 
with less numerical accuracy. In fact, we stop the time development of F of Eq. 
Q2TUJ1 before it reaches around 10 -3 in practical calculations. On the contrary, in 
the original ACCC method only the eigen-energy values in the bound state region of 
5 > 5o are used, and hence high numerical accuracy is needed for the extrapolation 
of the Pade rational function up to 5 = which is far from the bound state region. 
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§4. Conclusion 

We proposed a new method beyond the ACCC method in which in order to 
determine the coefficients of Pade approximant to the complex eigen energy we adopt 
the negative and positive energies which are obtained in bound state approximation 
or an extended version of bound state approximation to the description of resonances. 
The determination of the coefficients of Pade approximant was made by the use of 
the frictional cooling method which has been widely used in the works with AMD 
theory. The extended version of bound state approximation was already utilized in 
the previous work by present authors, and here we gave a full explanation of the new 
method. The new method enables us to extract the resonance component from the 
energy eigen-states in the bound state approximation which consist of the mixture 
of the resonance state and the continuum state components. 

We investigated the validity of the new method by the application to the second 
2+ state of 12 C which was discussed with the original ACCC method in our previous 
paper. We showed that the resonance states which were obtained in the extended 
version of bound state approximation correctly give the resonance energies which are 
determined by the use of the original ACCC method. The use of energy levels in 
the wide region of coupling constant which includes both resonance state and bound 
state was found to lead to the unique construction of the Pade rational function, 
which qualitatively well agrees with that of the original ACCC method. Thus we 
obtained similar values of energy and width to those given in the original ACCC 
method. The number of adopted eigen-energy states to construct the Pade rational 
function can be taken irrespectively of the dimension of the function, i.e. {N + M), 
while in the original ACCC method the number is fixed to (N + M). It should 
be noted that this new method needs much less numerical accuracy compared with 
the original one in which eigen value problem must be solved with extremely high 
numerical precision. We can conclude that our new method has several advantages 
over the original one and can be applied directly to resonance problems whenever 
the original ACCC method is applicable. 
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